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ABSTRACT 

SN 2006gy was the most luminous SN ever observed at the time of its discovery 
and the first of the newly defined class of super luminous supernovae (SLSNe). The 
extraordinary energetics of SN 2006gy and all SLSNe (> 10®^ erg) require either atypi- 
cally large explosion energies (e.g., pair-instability explosion) or the efficient conversion 
of kinetic into radiative energy (e.g., shock interaction). The mass-loss characteristics 
can therefore offer important clues regarding the progenitor system. For the case of SN 
2006gy, both a scattered and thermal light echo from circumstellar material (CSM) 
have been reported at later epochs (day ~ 800), ruling out the likelihood of a pair- 
instability event and leading to constraints on the characteristics of the CSM. Owing 
to the proximity of the SN to the bright host-galaxy nucleus, continued monitoring 
of the light echo has not been trivial, requiring the high resolution offered by the 
Hubble Space Teleseope (HST) or ground-based adaptive optics (AO). Here we report 
detections of SN 2006gy using HST and Keck AO at ~ 3000 days post-explosion and 
consider the emission mechanism for the very late-time light curve. While the optical 
light curve and optical spectral energy distribution are consistent with a continued 
scattered-light echo, a thermal echo is insufficient to power the AT'-band emission by 
day 3000. Instead, we present evidence for late-time infrared emission from dust that 
is radiatively heated by CSM interaction within an extremely dense dust shell, and 
we consider the implications on the CSM characteristics and progenitor system. 

Key words: circumstellar matter — supernovae: general — supernovae: individual 
(SN 2006gy) — dust, extinction — infrared: stars 


1 INTRODUCTION 


At the time of its discovery, supe rnova (SN) 2006kv was the 
most l uminous SN ever o bserved (lQuimbvl2006l : |Prieto et al.l 
I2OO6I : lOfek et ^ 20M_ Smith et al.l l2007l) .~The Type Iln 
spectrum (see iFilippenkd 19971 for a review of SN spectral 
classification) indicated that interaction with dense circum¬ 
stellar material (CSM) might play an important role in its 


energetics, but interpreting this event was not straightfor¬ 
ward. Aside from its peak luminosity, SN 2006gy has sev¬ 
eral other distinguishing characteristics. For example, the 
SN took nearly 70 days to rise to peak, remained brighter 
than —21 mag for 100 days, and had a total radiated 
energy of > 10®^ erg. 

A number of studies attempt to explain the early-time 
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Table 1. HST Observations 


UT Date 

Epoch 

Program 

PI 

Instrument 

Grating/Filter 

Central A 

Exp. 

Magnitude 


(days) 





(A) 

(s) 

(Vega) 

20121216 

2303 

13029 

Filippenko 

WFC3/UVIS 

WFC3/UVIS 

F625W 

F814W 

6241 

8026 

510 

680 

21.404 (0.013) 
20.948 (0.014) 





WFC3/UVIS 

F275W 

2707 

846 

> 25.3 

20130302 

2379 

13025 

Levan 

WFC3/UVIS 

F390W 

3922 

932 

22.882 (0.017) 





WFC3/IR 

F160W 

15369 

206 

> 14.7 





STIS/CCD 

G750L 

7751 

5228 

— 

20141114 

3001 

13287 

Fox 

STIS/CCD 

G430L 

4300 

5848 

— 

STIS/MAMA 

MIRNUV 

2305 

697 

> 27.0 





STIS/MAMA 

MIRFUV 

1452 

1957 

> 24.8 


properties of both SN 2006gy and the new class of “su- 
perluminous supernovae” (SLSNe) that emerged following 
its discovery, focusing on either a standard core-collapse 
event within dense CSM or the radioactive decay from sev¬ 
eral solar m asses of ^^Ni ge nerated in a pair-instability ex¬ 
plosion (see iGal-YamI |2012| . and references therein). Late¬ 
time observations offer an opportunity to rule out the 
pair-instability model with energetic constraints. Limited 
ground-based adaptive optics (AO) and Hubble Space Tele¬ 
scope (HST) photometry at ~ 300-800 days are signifi¬ 
cantly lower than predicted by the > 10 M© of ®®Ni pair- 
instability model th at would have been required to power th e 
light curve at peak llNomoto et al.|[2007l : ISmith et al.ll2008ll . 
While the day ~300 photometry is consistent with radioac¬ 
tive heating from a minimum of 2.5 Mq of ®®Ni, subsequent 
optical and infrared (IR) photometry reveals an observed 
excess of emission over the expected radioactive decay rate, 
sugges ting that the decl ine cannot be explained by ®®Co 
alone (iMiller et al.l I 2 OIOII . Given the seeming lack of evi¬ 
dence f or CSM interaction at earlier times, amongst other 
reasons. ISmith et al.l (l2008li propose the most likely scenario 
is a thermal-IR echo from a massive (~ 0.1 Mq) dust shell 
located ~ 1 ly from the SN that was radiatively heated by 
the SN peak luminosity, but they cannot rule out a large 
®®Co luminosity that was reabsorbed by dust. iMiller et al.l 
1 I 2 OIOII later eliminate a large ®®Co mass with Keck AO ob¬ 
servations that show a slower decline in the IR luminosity 
than expected from radioactive decay. A relatively blue late¬ 
time optical colour from those epochs also suggests the pres¬ 
ence of a scattered-light echo. These results rule out various 
pair-instability SN models. 

Thermal li ght echo mo dels predict the evolution of the 
IR light curve IIDweklflQSSll . A scattered-light echo should 
have an optical spectrum reminiscent of the SN spectrum 
at peak (as opposed to late-time CSM interaction). At only 
~ 1 . 2 " from the bright galaxy nuclear bulge with a promi¬ 
nent dust lane, however, late-time (> 300 day) data for SN 
2006gy have been relatively difficult to obtain. No spectrum 
of SN 2006gy has been published since day 237. (Smith et al.l 
ll2008ll published a two-dimensional (2D) ground-based Keck 
spectrum of SN 2006gy on day 364, but a reliable ID ex¬ 
traction was too difficult given the overwhelming galactic 
nucleus.] 

Here we present ongoing observations of SN 2006gy 
through day 3024 to test for the presence and expected 
evolution of both a scattered and thermal light echo. Data 


include day 3024 Keck/NIRC2-AO-LGS A'-band photom¬ 
etry, a day 3001 HST/STIS spectrum of SN 2006gy, and 
HST/WFCS and STIS visible and ultraviolet (UV) photom¬ 
etry. Section [ 2 ] presents the observations. Sections [3] and [4] 
analyse the data in the context of a light echo, while Section 
[S] considers the case of CSM interaction. Section presents 
the implications of the results and summarises our conclu¬ 
sions. 

Throughout this paper we assume that the distance to 
NGC 1260 (the host galaxy of SN 2006gy) i s 73.1 Mpc (red- 
shift a = 0.018), and following ISmith et ahl ll2007l') we adopt 
E{B — V) = 0.54 mag as the reddening toward SN 2006gy 
within the host galaxy, while Galactic extinction accounts 
for E{B — V) = 0.18 mag, leading to a total colour excess 
toward SN 2006gy of E{B — V) = 0.72 mag. Unless oth¬ 
erwise noted, all spectral energy distributions (SEDs) and 
spectra have been corrected for this colour excess assum- 
ing Ry = Ay / E{B — U) = 3.1 using the reddening law of 
ICardelli et al.l lll989h . 


2 OBSERVATIONS 

2.1 Hubble Space Telescope 

Table [1] lists details regarding new data on SN 2006gy ob¬ 
served by HST as part of programs GO-13287 (PI Fox), 
GO-13025 (PI Levan), and GO-13029 (PI Filippenko). 


2.1.1 WFC3 Photometry 

SN 2006gy was observed with the HST/WFCi UVIS and 
IR channels, summarised in Table [T] and displayed in Fig- 
ure[T] Photometry was extracted from the individu al WFC3 
“Jlt. fits” images in all bands using Dolphot v2.0 (IPolphinl 
I 2 OOOII . The input parameters are those recommended by the 
Dolphot WFC3 Users’ Manual. Aperture corrections were 
applied. The region in the F160W filter is too confused for 
accurate photometry. Instead, the upper limit is based on 
unreal sources that Dolphot thinks it has detected in the 
main body of the galaxy, near the nucleus and the dust 
lane. Table [T] lists the resulting magnitudes in the WFC3 
flight system (Vegamag). 
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HST/WFC3/F390W 

20130302 


HST/WFC3/F160W 

20130302 


HST/WFC3/F814W 

20121216 


SN 2006gy 


HST/WFC3/F625W 

20121216 


Figure 1. HST imaging of SN 2006gy. Details of the instrument, filter, and observation date are labeled in each image and listed in 
Table [T] 


2.1.2 STIS Photometry 

SN 2006gy was observed with the HST/STIS NUV and FUV 
MAMA channels, summarised in Table [T] and displayed in 
Figure [T] No source is detected at the position of the SN (a 
bright knot in one the galactic arms is all that is visible). We 
calculate upper limits using standard photometry techniques 
in IRAF’s APPHOT package and convert to magnitudes us¬ 
ing procedures outlined by the STIS Data Handbook Section 
5.3. 

2.1.3 STIS Spectroscopy 

Table m summarises the optical spectra of SN 2006gy ob¬ 
tained with HST/STIS. Given that the SN magnitude was 
unknown at the time of the observations, the SN was cen¬ 
tred in the 52 x 0.2” slit by offsetting from a nearby star. 
The location of the supernova in the 2D image is identified 
by matching a simulated cross-dispersion (XD) profile from 
WFC3/F814 observations (day 2303) with the actual XD 


profile of the STIS/G750L spectrogram. We choose these 
two observations for comparison since they are most closely 
matched in wavelength coverage and include the prominent 
Ha line. 

Although we specified the slit position on the sky in 
the original Phase 2 Hies, we consider a range of actual 
slit positions. For each position, the WFC3 XD profile is 
simulated by summing the WFG3 pipeline calibrated image 
(_drz.fits) in the dispersion direction over 5 pixels, which 
corresponds to the STIS slit width on the WFC3 detector. 
The STIS XD profile is created by collapsing the 2D spec¬ 
trogram in position space and rebinning to the WFG3 plate 
scale. The slit position is then constrained by matching the 
WFG3 and STIS XD profiles (see Figure [D). 

The ID spectrum for each observation is extracted us¬ 
ing the GALSTIS custom extraction software stistools.x lc0. 


^ http://ssb.stsci.edu/doc/stsci_python_x/stistools.doc/html/xld.html7highlight 
stistools.xld 
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Normalized XD profiles from STIS (NUV) and WFC3 



Figure 2. Slit position orientation confirmation, (left) The ifS'T/WFC3/F814W image with the planned slit orientation highlighted, 
(right) A simulated cross-dispersion profile (blue) created by summing rows within the slit on the left-hand side. The green line 
represents the actual STIS cross-dispersion profile, calculated by median-stacking the columns in the pipeline-reduced STIS spectrum 
(ocdd04010_crj.fits). The similarity of these two profiles indicates the SN, highlighted by the grey bar just south of the galaxy nucleus, 
is centred in the slit. 


The default extraction parameters for STIS are defined for 
an isolated point source. For both the G430L and G750L the 
default extraction box width is 7 pixels and the background 
extraction box width is 5 pixels. SN 2006gy, however, is quite 
faint and located just a few pixels in position space below 
the galaxy (0.05”pixel“^; Figure^). We therefore reduce the 
extraction box width to 3 pixels to optimise the signal-to- 
noise ratio (S/N). After locating the SN position along the 
slit above, the optimal extraction pixel position was deter¬ 
mined by shifting the extraction box along the 2D spec¬ 
trogram in position space and optimising the Ha S/N. The 
three extractions above and below this SN position were me¬ 
dian combined to estimate the background. All of the STIS 
spectra are combined to produce a single spectrum using 
the splice tool in STSDA^. We ignore the edge columns by 
setting their data-quality flags to 4. Figured plots the final, 
background-subtracted spectrum. 


present in all the images was used to align and, ultimately, 
coadd the dithered images. 

The resulting data were analysed using standard IR 
analysis techniques utilising SExtractoiQ. Given the steep 
gradient of the underlying galax y, we chose to subtract this 
contribution by using GALFIT (IPeng et al.ll2003 . 120101 ) to 
model NGG 1260 with a radial Sersic 1 19631 ) profile. We 
use a bright, isolated star as a model for the point-spread 
function. We allowed the exponent n, half-light radius, axis 
ratio, position angle, galaxy position, and sky background to 
vary. Removing the principal source of variation in the back¬ 
ground yields a more robust local background subtraction. 
Galibration was p erformed using field s tars with reported 
fluxes in 2MASS (ISkrutskie et al.ll2006l). Table [2] l ists the 
new photometry and that reported bv iMiller et al.l ll2010l) . 
Figure [6] plots this photometry along with the WFC3 pho¬ 
tometry presented in Table [T] 


2.2 Keck/NIRC2-AO-LGS 

SN 2006gy was observed with Keck/NIRC2-AO-LGS 
dWizinowich et al]l2006l ) on 2014 Dec. 7 UT using the wide- 
held camera and K' filter. The complete set of observations 
include two repeats of the “bxy9” pattern with 1" dithers, 
where each frame consisted of three coadded 8 s exposures 
(i.e., 24s frame). To reduce the data, we subtracted a me¬ 
dian bias frame, applied fiat-field corrections to each expo¬ 
sure, and corrected for astrometric distortiorf^. A bright star 

^ http://stsdas.stsci.edu/cgi-bin/gethelp.cgi?splice 
® https://www2.keck.hawaii.edu/inst/nirc2/dewarp.html 


2.3 Keck Optical Spectroscopy 

For a comparative analysis in Section[5l we also present pre¬ 
viously unpublished optical spectra of SNe Iln 2005ip and 
2010jl at similar late-time epochs, summarised in Table [S] 
The spectra were obtained with t he dual-arm Lo w Resolu¬ 
tion Imaging Spectrometer (LRIS: [Oke et al1ll995l) mounted 
on the 10-m Keck I telescope with the slit aligned along 
the parallactic a ngle to minimise differential light losses 
llFilippenkdl 19821) . The spectra were reduced using standard 

^ SExtractor can be accessed from 

http: //www.astromatic.net/software 
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Figure 3. The 2D spectrum of SN 2006gy obtained with the G750L grating (ocdd04010_crj.fits). The SN is highlighted by the faint 
trace just below the bright galaxy nucleus. The Ha line, highlighted by the inset, is the most obvious feature. 



Rest Wavelength (A) 



5000 6000 7000 8000 9000 10000 

Rest Wavelength (A) 

Figure 4. Background-subtracted i^S'T/STIS spectrum of SN 2006gy on day 3001 post-explosion, corrected for extinction assuming a 
reddening of E{B — V) = 0.72 mag. The light blue plots the unsmoothed spectrum, while the dark blue plots the spectrum that has 
been 10-pixel boxcar smoothed. Grey bars signify pixels flagged for having bad data quality. The spectrum is relatively featureless, flat, 
and noisy. We detect Ha (shaded red), but no other obvious lines can be identified. 
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Table 2. Keck AO Observations of SN 2006gy 


Date 

(UT) 

Epoch 

(days) 

Filter 

Mag 

(Vega) 

2007 Sep. 29 

398 

A' 

14.91 ± 0.17 

2007 Dec. 2 

461 

H 

16.8 ± 0.3 

2007 Dec. 2 

461 

K' 

15.02 ± 0.17 

2008 Aug. 25 

723 

K' 

15.59 ± 0.21 

2014 Dec. 07 

3024 

K' 

18.10 ± 0.17 


Table 3. Summary of Keck/LRIS Optical Spectra 


SN 

JD - 

Epoch 

Res. 

Exp. 


2,450,000 

(days) 

(A) 

(s) 

2005ip 

6778 

3024 

~ 6 

1200 

2010jl 

6778 

1290 

~ 6 

1200 


techniques (e.g., iFolev et alJ 1200^ : ISilverman et al.|l2012l) . 
Routine CCD processing and spectrum extraction were com¬ 
pleted with IRAI0, a nd the data were extracted with the 
optimal algorithm of iHorn 3 lIlQSdl . We obtained the wave¬ 
length scale from low-order polynomial hts to calibration- 
lamp spectra. Small wavelength shifts were then applied to 
the data after cross-correlating a template-sky spectrum to 
an extracted night-sky spectrum. Using our own IDL rou¬ 
tines, we fit a spectrophotometric standard-star spectrum to 
the data in order to fl ux calibrate the SN and to remove tel¬ 
luric absorption lines (IWade fc Hornelll988l : lMatheson et al.l 

l200Cil . 


3 A SCATTERED OPTICAL LIGHT ECHO? 

The observed late-time optical emission from SN 2006gy 
has been attribut e d previously to a scattered-light echo 
llSmith et al.l l2008l : iMiller et akl l201Cll l , which is a product 
of scattered light from the SN light curve emerging from a 
paraboloid of revolution with the supe rnova as its focus and 
its axis along the line of si ght fe.g.. I^de fc Evans! Il9^ : 
lDweklll983l : IChevalierlfloS^) . We consider this scenario in 
the context of the new data presented in this article. 

Figure [7] plots the optical SED of SN 2006gy on day 
2379, which corresponds to the single epoch of F390W ob¬ 
servations. Although data through the other filters were not 
obtained at this particular epoch, we extrapolate the F625W 
and F814W photometry shown in Figure [S] This figure goes 
on to compare the photometry to the expected scattered 
echo spectrum, which is the cumulative scattering of the 
whole light curve. To simulate this spectrum, we construct 
a synthetic integrated spectrum by mean-combining the in¬ 
dividual spectra from day 36 (pre-peak), day 71 (peak), day 
122, and day 177 (spaced roughly every 40 davs: ISmith et aP 
|2010H . We then assume a A wavelength dependence 

® IRAF: the Image Reduction and Analysis Facility is distributed 
by the National Optical Astronomy Observatory, which is oper¬ 
ated by the Association of Universities for Research in Astronomy 
(AURA), Inc., under cooperative agreement with the US National 
Science Foundation (NSF). 



Figure 5. Keck/NIRC2/LGS-AO image of SN 2006gy obtained 
on 2014 Dec. 7 UT. 

for the scattering , which can be considered typical (e.g., 
iMiller et al.ll201(l ). Since the flux of the scattered-light echo 
depends on the specihc arrangement of the dust, the pre¬ 
cise scale factor is unknown. Instead, Figure [7] scales the 
spectra to the F814W photometry. Overall, the synthetic 
photometry for a scattered-light echo with a wave¬ 

length dependence is consistent with the observed SED on 
day 2379. 

4 A THERMAL-IR LIGHT ECHO? 

4.1 The Peak of the Thermal SED 

Little late-time IR colour information exists for SN 2006gy 
given that the observations require high-resolution ground- 
based AO or HST photometry. Prior to this article, just 
1 H- and 4 A'-band observations existed at > 1 yr post¬ 
explosion. The equilibrium dust temperature and, thereby, 
the peak of the corresponding SED rem ain relatively uncon¬ 
strained. Both ISmith et al.l |2008h and IMiller et al.1 (|2nifll ) 
place a lower limit on the dust IR luminosity by assuming 
that (1) all of the near-IR luminosity can be attributed to 
a thermal light echo, and (2) the near-IR luminosity peaks 
in the K' band given the very red H — K' colour observed 
on day ~ 400. Here we present a quantitative analysis of 
the A'-band contribution to the total IR luminosity in the 
thermal-IR echo scenario. 

First, we calculate the equilibrium temperature of dust 
over a range of distances f rom the S N fol lowing dust 
heating models outlined by IFox et ;il] (120101) . Assuming 
0.1 qm graphite dust grains, Figure[8Ka) plots the dust tem¬ 
perature as a function of the central energy source (i.e., SN 
peak luminosity) for dust at distances ranging from 0.2 to 
4 ly. For any given luminosity, the dust is heated to ex¬ 
pectedly lower temperatures at larger distances. Figure ISTb) 
plots the dust temperature as a function of distance for the 
specihc case of the observed peak luminosity of SN 2006gy 
(Apeak ~ Lq). 
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SN 2006gy Photometry 



Epoch(days) 

Figure 6. Light curve of SN 2006gy through day 3024, including data from ISmith et all ll2nn8h and iMiller et alJ ll2ni(]l i. Dotted lines 
for the F625W and F814W filters beyond day ~ 2300 illustrate the extrapolations used to calculate the photometry plotted in Figure[7l 
corresponding to the single epoch of F390W photometry on day 2379. Offsets are applied only for plot clarity. 



3000 4000 5000 6000 7000 8000 9000 10000 

Rest Wavelength (A) 

Figure 7. The optical SED of SN 2006gy on day 2379, which corresponds to the single epoch of F390W observations. Overplotted in 
grey (solid) is the synthetic integrated spectrum of SN 2006gy around peak light, constructed by mean-combining the day 36 (pre-peak), 
day 71 (peak), day 122, and day 177 data. A wavelength dependence is assumed for the scattered-light peak spectrum (dashed 

grey). Both spectra are scaled to the F814W photometry. The synthetic photometry of the scattered spectrum is consistent with the 
observed SED. 


Since the dust temperature can be written as a function 
of radius from the SN, we can also plot the fraction of the 
total IR luminosity emitted in K' as a function of radius 
(Figure [HI). At larger distances, where the dust temperature 
drops and the SED peak shifts to longer wavelengths, the 
fraction of the total IR flux emitted in K' decreases. 

For the case of an IR echo, the quantitative relation¬ 
ship between observation epoch and the emitting dust shell 
radius is not straightforward because the paraboloid inter¬ 
sects sm^ fractions of many thin shells at any given instant 
(see IDw^I 1983I ). While modeling the integrated flux from 
the many thin shells is beyond the scope of this paper, the 
hottest dust at any epoch, t, is located at a radius R = ct/2, 
where c is the speed of light. This radius therefore sets the 
upper limit to the fractional K'-hand emission (all other 
contributing shells in the paraboloid have lower tempera¬ 
tures). By writing the observation epoch as a function of 
the hottest dust shell radius, t = 2R/c, Figure [9] also plots 
the maximum fraction of the total IR luminosity emitted 


in K' as a function of time post-explosion. The fraction of 
K'-hand flux at a given epoch is independent of the CSM 
density or geometry. At early times, the K' flux represents 
only ~ 15% of the total IR flux, and by late times, this 
fraction drops to <1%. 

4.2 Other Potential W'-Band Flux Sources 

Besides the equilibrium thermal-IR emission, other potential 
emission sources may contribute to the Tf'-band flux at late 
times, including (1) the scattered-light echo at 2 pm, (2) 
thermal emission from hotter, smaller grains that are not in 
thermal equilibrium, and (3) H 2 line and CO band emission. 
Here we consider possible contributions from these sources. 

Scattered-Light Echo at 2 pm: Figure (3 approxi¬ 
mates the spectrum and SED of a scattered-light echo. Fig¬ 
ure [TU] goes on to extend this SED into the IR by fitting the 
SED with a blackbody and extending the Rayleigh-Jeans tail 
to the K' band. The expected fraction of the K'- to R-band 
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Figure 8. Relationship between dust mass, dust temperature, and distance from the heating source, (left) The dust temperature as 
a function of the central energy source (i.e., SN peak luminosity) for 0.1 \im graphite dust grains at distances ranging from 0.2 to 4 ly 
(dotted lines). Overplotted as a reference are the luminosities corresponding to both the peak and late-time plateau of SN 2006gy. (right) 
The dust temperature as a function of distance for the specific case of the observed peak luminosity of SN 2006gy (Tpeak ~ 10^^ Lq). 



Figure 9. The fraction of K' to total IR luminosity as a function 
of dust temperature, which can be written in terms of radius 
(Figure or time (t = 2R/c), assuming a thermal light echo 
model and 0.1 pm grains. 


Figure 10. The scattered optical light echo from Figure [7| ex¬ 
tended into the IR assuming a Rayleigh-Jeans tail. For a light- 
echo scenario, the i^'-band contribution would be insignificant 
compared to what is observed. 


flux from an unreddened scattered-light echo would be only 
0.02. From Tables^and[^ we calculate the observed ratio 
on day 2303 (by interpolating the K'-hsiiid fluxes and assum¬ 
ing a linear decline): K'jR = 17.3 mag/20.9 mag = (9.2 x 
10“^® erg cm“^ A"^)/(6.7 x 10"^® erg s"^ cm"^ A~^) « 
1.38. We therefore rule out any significant contribution from 
the light echo at 2 (rm in our models. 

Small, Hot Grains: Small grains do not radiate as 
blackbodies. At a given distance from the SN, smaller grains 
will therefore be hotter than the equ ilibrium blackbody tem - 
perature exhibited by larger grains. iTemim &: DwekI ll2013l l 
show in their Figure 4, however, that the dust grain temper¬ 
ature peaks and plateaus for dust grain sizes a < 0.1 pm (as¬ 
suming a constant heating source). Our models in Section 
14.11 already assume grain sizes a — 0.1 pm and include the 
as sociated absorp tion and emission coefficients (see details 
in iFox et al.ll201(]l l. The modeled fraction of the IR luminos¬ 


ity emitted in is''-band flux (Figure E} therefore represents 
only an upper limit. 

H 2 line and CO band emission: Finally, we note 
that additional K' emission may originate from either 
H 2 line or CO band emission, but we assume the contribu¬ 
tion is negligible in a broadband filter and do not consider 
these contributions in our models. 

4.3 Energy Budget 

The integrated optical energy output from the SN photo- 
sphere throughout the first ~ 200 days is ~ 2.5 x 10®^ erg 
JSmith et aLlboiol l. Assuming that the IR luminosity peaks 
in the K' band and a constan t K' luminosity of 2 x 10® L© for 
600 days, iMiller et al.l ll2010ll estimate the total emitted IR 
energy at Eir > 4 x 10"*® erg. Making a similar assumption 
about the thermal emission peak wavelength and integrating 
over the observed K' light curve in Figure El we calculate a 
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Figure 11. For the light-echo model, the fraction of K'- to R- 
band flux as a function of dust temperature, which can be written 
in terms of radius (Figure (8)1 or time (t = 2iJ/c) Overplotted in 
red are the actual observed value from Figure |6] 


similar total emitted energy through day 3000 (it turns out 
the assumption of a Jf'-band luminosity of 2 x 10® Lq for 
600 days was an overestimate). 

Figure [5] shows, however, that in the light-echo model 
the fraction of the IR luminosity emitted in the K' band 
is 15% at day 400 and < 1% by day 3000. Accounting for 
the fractional output in K', the lower limit on the total IR 
luminosity is actually Eib. ^ 4 x 10®^ erg. This calculation 
sets only a lower limit because the fractional A'-band output 
represents only the hottest dust shell. The total radiated IR 
energy from the putative echo is therefore greater than the 
total SN output, which is not even possible for a case of an 
optically thick shell. This energetics argument alone suggests 
that the IR echo argument is unfeasible at day 3000. 


4.4 R-to-K' Band Ratio 


We also consider the colour evolution in the context of the 
light-echo scenario. Specifically, we derive the ratio of the R 
and K' bands, which are both observables. The scattered op¬ 
tical and thermal-IR fluxes can be approximated as a func¬ 
tion of radius: 


^phot Atphotrscat/" (^) 

2R/c ’ 


Lr{R) = 
and 


Y / _ Aphot Afphot'rd tjRi 

TrTc 


( 1 ) 

( 2 ) 


where Lphot is the photosphere luminosity over a given time 
Afphot, Tscat is the dust scattering coefhcient, ra is the 
dust absorption coefficient, f{6) is the fraction of forward- 
scattered electrons, and Lr/ /Lir is the calculated IR frac¬ 
tional output in K' (see Figure O- The ratio of the R- to 
A'-band flux versus radius can then be written as 



Tscat/(^) ^ Lir 

Td Tk' 


C X 


Lir 

Lk' 


(3) 


While the values for Tscat, Td, and f{6) may require detailed 
derivations, they can all be considered constants for the pur¬ 
poses of this analysis, since the ratio of these values will not 
change as a function of the dust-shell radius. 

Figure fTT] combines Equation [3] with the analysis from 
Section r4.3l to plot the fraction of the R- to A'-band flux as a 
function of radius. As the echo shifts to shells at larger radii, 
the fraction of thermal echo emitted in K' declines while 
the scattered optical light echo remains constant. Overplot¬ 
ted are the measured ratio of the R- and A'-band fluxes 
from Figure [HI For the first epoch, the ratio may be consis¬ 
tent with the light-echo scenario (depending on the values 
of Tscat, Td, and f{0)). The trend over ~ 1500 days, how¬ 
ever, does not follow the predicted shift toward a signifi¬ 
cantly larger optical flux. These calculations again suggest 
that the thermal-IR echo scenario is unlikely to power the 
light-curve plateau, particularly by day 3000. 


4.5 Dust Temperature and Mass 

iMiller et al.l ll2010li place a lower limit on the dust shell in¬ 
ner radius, Ri > 1.5 x 10^® cm, based on both the beginning 
of the near-IR echo plateau and the SN peak luminosity. 
This radius should be reconsidered for two reasons. First, a 
distance of Ri > 1.5 x 10^® cm implies a light-travel time 
of > 1.5 yr. This is inconsistent with their reported near- 
IR excess from a thermal echo as e arly as day 130. Sec ond, 
for the peak luminosity argument, iMiller et al.l (l2010li as¬ 
sumed a dust va porisat i on te mperature of T-u « 1000 K, 
similar to that of lPwekl lll983ll . In reality, the dust vapori¬ 
sation temperature for sili con and graphite dust particles is 
1500 ^ ^ 2000 K (e.g.. [Gall et al.ll2014f). Given t he dust 

shell radius dependence in Equation 2 of Miller et al. IH^, 
Ri oc the calculated vaporisation radius decreases by 

a factor of 3-6, or Ri > (2.5-5) x 10^^ cm. 

This new estimate of Ri changes some of the analy¬ 
sis and interpretation of SN 2006gy. First, when taking into 
account the light-travel time, this radius is now consistent 
with an observed near-IR excess as early as day 130. Sec¬ 
ond, Equation 3 of iMiller et al.l (l201(]l l now yields a different 
estimate of the total dust mass in the IR-echo scenario, as¬ 
suming a r~^ wind. If the outer dust radius is given by the 
light-travel time from the most recent set of Keck/NIRC2 
observations (day 3024), then i ?2 ^ ct/2 = 4 x 10*^® cm = 
(8-16)7?i. According to Equation 3 in IMiller et al.l (l2010l) . 
this increases the estimate of the total dust mass by nearly 
a factor of ~ 2-3, or Md « 0.2-0.3 Mq. 

This dust mass derived from the light-echo scenario is 
difficult to reconcile with t he observed A'-ba nd m agnitudes. 
Follow ing the analysis of IFox et ahl (l2010ll and IFox et ahl 
ll201lD , a single A'-band flux can be fit as a function of 
the dust mass and temperature (e.g., a higher tempera¬ 
ture requires a lower dust mass and vice versa). Figure \TI\ 
shows this relationship for the measured A'-band flux at 
each epoch. 

On day 3024, Figure [5] indicates that for the light-echo 
model, the inner dust-shell radius, R ~ 5x 10^® cm, will have 
have a temperature of <300 K. The hottest dust serves as a 
useful lower limit on the required dust mass because colder 
dust would contribute even less flux to the measured A'- 
band flux. Figure [12] reveals that to recreate the observed 
A' flux with 300 K dust would require > 10^ M©. This mass 
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Figure 12. A single iT'-band flux can be fit as a function of 
the dust mass and temperature (a higher temperature requires a 
lower dust mass and vice versa). The different combinations from 
this relationship are calculated for the measured K' flux at each 
epoch. 


is significantly greater than the dust mass derived above 
from the light-echo equations. Furthermore, such a large 
dust mass would require > 10"* Mq of gas (assuming typical 
gas-to-dust mass ratios of ~ 100). Taken all together, these 
results offer further evidence that the thermal-echo scenario 
is not likely the dominant mechanism powering the late-time 
IR plateau of SN 2006gy by day 3000. 


5 IR EMISSION FROM GSM INTERACTION? 

Other physical scenarios may explain the late-time IR 
plateau of SN 2006gy, including radioactive decay from a 
pair-instability SN, new dust formation, collisi onal heating, 
and r adiat ive heating f rom C SM interaction. ISmith et al.l 
(l2008l) and iMiller et al.l (|2010| ) rule out most of these scenar¬ 
ios with earlier data. Specifically, ISmith et al.l (l2008l ) argue 
against CSM interaction at <800 days given the absence of 
broad or intermediate Ho emission in a day 364 2D spec¬ 
trum (but see Section [5.51 below), weak X-rays, and no radio 
detection. They also rule out the possibility of dust obscu¬ 
ration, since the calculated dust mass would be insufficient 
to hide the expected Ha luminosity. These calculations as¬ 
sume a dust temperature of Td = 1300 K so that the energy 
distribution peaks in the K' band. 

While the light echo may have dominated at these ear¬ 
lier epochs, such assumptions are no longer valid by day 
3000. Here we consider the possibility that CSM interac¬ 
tion dominates the IR light curve at these later epochs with 
new assumptions and data. In this scenario, a pre-existing 
dust shell is radiatively heated by X-ray and /or optical emis¬ 
sion g enerated by ongoing CSM interaction llFox et al.ll201ll . 
l2013h . Throughout the analysis, we assume a spherically 
symmetric distribution of dust. 


5.1 Dust Temperature and Mass 

Unlike the thermal-echo scenario, where only a fraction of 
each dust shell emits at a given epoch, the observed IR flux 


in the CSM interaction scenario corresponds to the entire 
dust shell (as long as the light crossing time is shorter than 
the duration of CSM interaction). Assuming spherical sym¬ 
metry, the IR flux can therefore be converted into a min¬ 
imum dust shell blackbody radius, rbb = (I/ir/Itto-T,!)^/^. 
Figure [12] plots the blackbody radius corresponding to each 
combination of dust mass and temperature used to fit the 
K' fluxes at each epoch (see Section r4.5ll . although a number 
of degeneracies remain. 

Independent constraints can break these degeneracies 
in Figure 1121 Figure [8] demonstrates that for a peak lu¬ 
minosity of 10^^ Lq, the dust vaporisation radius (as¬ 
suming a vaporisation temperature of Td « 1500 K) is 
~ 2 X 10^^ cm = 0.067 pc. For rbb = 0.067 pc. Figure 
13 shows a corresponding dust-shell temperature and mass 
on day 3024 of ~ 450 K and 8 x 10“^ Mq, respectively. 
Figures [8] and [9] highlight that a temperature Td « 450 K 
corresponds to a fractional A'-band output of 2%. (These 
calculations assume a dust-shell radius equal to the vapor¬ 
isation radius, which requires a constant pre-SN mass loss. 
Of course, the pre-SN mass loss may not have been con¬ 
stant and the dust-shell radius may be larger. We consider 
this possibility later.) 


5.2 Energy Budget 

Accounting for this fractional IR flux in K', the total emit¬ 
ted IR and optical energy throughout the plateau phase is 
Air > 2 X 10®° erg and Aopt ~ 2 x 10'*° erg (see Section IT3| . 
For the CSM interaction scenario, the thermal-IR emission 
results when the dust shell absorbs and reradiates the opti¬ 
cal flux generated by the shocks. The implied optical depth 
is therefore 


Td 


— 1 X In 
2.3. 


E. 


opt 


Eopt + Air 


(4) 


Assuming most of the optical energy is emitted in the 
R band, the extinction can be estimated as An k, Td = 
2.3, w hich is slightly larger than estimated bv ISmith et al.l 
(l2008l) . A slightly warmer dust shell or additional visible 
radiation outside of R would decrease the value of An toward 
the value of 1.5 measured bv ISmith et al.l (|2008l ). 

The integrated optical energy output from the SN pho¬ 
tosphere throughout the first ~ 200 days is ~ 2.5 x 10®^ erg 
(see Section 14.311 . Assuming CSM interaction dominates at 
later epochs, the total radiated energy emitted from the SN 
throughout day 3000 is therefore Atot ~ 2.5 x 10®*^ -f 4 x 
10®° -I- 2 X 10"^° = 2.9 X 10®^ erg. This number is larger 
than the radiated energy output by most SNe, but consis¬ 
tent with many SLSN, especially considering observations of 
SN 2006gy ex tend all the way through day 3000 (see Table 
1 in iGal-YamllioIi) . 


5.3 R-to-K' Band Ratio 

Figure [TT] plots the ratio of the observed R- and A'-band 
fluxes. The ratio remains nearly constant throughout the 
extent of the observations. While this trend is not consis¬ 
tent with the expected colour evolution of a light echo, it is 
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Figure 13. The optical SED of SN 2006gy on day 2379, which corresponds to the only epoch of F390W observations (similar to Figure 
[ 7 ) 1 . Overplotted is the day 1290 spectrum of SN 2010jl (red) and the day 3024 spectrum of SN 2005ip (red), both which show evidence 
for CSM interaction. Both spectra are scaled to the F814W photometry. Strong Ha for both SNe 2005ip and 2010jl results in a poor 
match with the observed F625W emission from SN 2006gy. 


expected for the CSM interaction scenario. If the shock re¬ 
mains interior to a relatively static dust shell and the optical 
depth stays constant, the ratio of the absorbed to emitted 
flux will also be constant. The K' flux may be expected 
to drop slightly as dust is destroyed by the forward shock, 
leaving only dust at larger radii. 


5.4 Can CSM Interaction Also Explain the 
Optical Emission? 

Figure [ 7 ] already shows that the observed SED of SN 2006gy 
on day 2379 is consistent with a scattered-light echo. We 
now consider the possibility that the optical SED may also 
be consistent with CSM interaction. Figure [13] again plots 
the observed SED of SN 2006gy on day 2379. This time, 
however, the figure goes on to compare the photometry to 
the late-time Keck spectra obtained of both SNe 2005ip (day 
3024) and 2010jl (day 1290), two Type Iln explosions that 
are known to exhibit significant CSM interaction at very late 
epochs. To make for a fair comparison, both spectra are red¬ 
dened by a similar amount as SN 2006gy {E{B — V) — 0.72 
mag). Since we expect minimal contributions of CSM inter¬ 
action at redder wavelengths, the spectra are scaled so that 
their synthetic photometry matches the F814W photometry 
of SN 2006gy. 

The synthetic photometry of the two spectra is compa¬ 
rable with the F390W and F814W photometry of SN 2006gy 
(SN 2010jl more so than SN 2005ip), but it significantly over¬ 
estimates the F625W photometry. The most likely explana¬ 
tion for this is that both SNe 2010jl and 2005ip have strong 
Ha emission lines that are consistent with stro ng and on¬ 
going CSM interaction (ISmith et al.ll200^ . I2OI2I) . while SN 
2006gy does not exhibit prominent Ha. We search more care¬ 
fully for Ha in the day 3001 STIS spectrum of SN 2006gy 
in Section ESI below. 

Since the SED and spectral shapes are at least some¬ 
what consistent, however, this comparison raises the ques¬ 
tion of whether some other late-time spectrum with less 
CSM interactio n may offer a better fit (see, for example, 
IFox et ^l2013l ). While this paper does not compare an ex¬ 


haustive set of late-time spectra exhibiting CSM interaction, 
we point out that it would be difficult (if not unphysical) 
to produce significant line emission in the F390W band¬ 
pass (e.g., “blue pseudo-continuum”) without also generat¬ 
ing strong Ha emission in the F625W filter. Furthermore, 
these late-time spectra may be misleading because the SNe 
can be in a star cluster, have an optical echo of their own, 
or have bad subtraction of galaxy light, which is common 
when the SNe are faint. We therefore conclude that the op¬ 
tical emission from SN 2006gy is most likely dominated by 
a scattered-light echo, although CSM interaction may con¬ 
tribute a small fraction of the observed emission. 


5.5 STIS Spectrum and the Ha Flux 

Figure |3] plots the STIS spectrum corrected for extinction 
assuming a reddening E{B — V) — 0.72 mag. While the day 
3001 STIS spectrum is noisy, we can identify the Ha line, 
particularly in the 2D image (see Figure|3]|. Figure lTTl further 
confirms the presence of a broad, full width at half-maximum 
intensity (FWHM) 2000±200 km s“^ profile. This velocity is 
somewhat co nsistent with the sp eed of the post-shock gas (~ 
4000 km s~^; [Smith et al.ll2008l) . which shouldn’t decelerate 
significantly if the shock is plowing into a r~^ wind. No other 
lines can be definitively identified. 

We estimate the total output from the Ha line by in¬ 
tegrating over the Ha line profile after subtracting a con¬ 
tinuum of ~ 5 X 10“^® erg s“^ cm“^ A“^. The inte¬ 
grated flux of the line before correcting for extinction is 
(6.5 ±0.7) X 10®® erg s“® which is about an or der of magni¬ 
tude l arger than the upper limit placed on it bv ISmith et al.l 
ll2008l) . Some of this emission, however, originates in the 
scattered-light echo (see Section 0. We approximate the 
contribution from CSM interaction by subtracting off the ex¬ 
cess SN flux in the F625W filter compared to the synthetic 
photometry of the scattered-light-echo model in Figure [T] 
which comes out to ~ 4 x 10®® erg s“^, or 5% of the F625W 
flux (which is also within the error bar). We point out that 
the STIS spectrum was obtained nearly two years after the 
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SED constructed in Figure[71 so the relative contributions of 
the scattered echo and CSM interaction may have changed. 

Despite the identification of the broad Hq line, we con¬ 
sider the possibility that there is also some narro w-line con¬ 
tribu tion from an underlying H II region (e.g., IFox et ^ 
l2013t) . For a comparison, we consider the H II region of the 
Carina Nebula, the host to t] Carinae. While large when 
considering our local Milky Way neighbourhood, the Neb¬ 
ula is small relative to other known massive star form¬ 
ing complexes, consistent with the home environment of 
the likely massive progenitor to SN 2006gy. At the dis¬ 
tance of SN 2006gy, the whole Nebula would span only 
0 .2", which would be spatially unresolved in these spec¬ 
tra. The Carina Nebula has an integrated Ha luminosity of 
10® Lq JSmith fc Brook^ l200if l. which would leave nearly 
none of the remaining available Ha line flux to originate 
from CSM interacti on. This flux is with in th e limits indepen¬ 
dently calculated bv ISmith et al.l (l2010li and lAgnoletto et al.l 
||2009() . 


5.6 Powering the IR Emission 


Although we do not have a constraint on the dust mass or 
shell radius. Figure [12] shows that the observed K'-h&nd flux 
on day 3024 can be generated with 450 K dust in a spheri¬ 
cally s^unmetoc shell at the vaporisation radius. Figure 8(b) 
from iFox et al.l (|2010| ) calculates that 0.1 pm graphite dust 
at this temperature and radius requires a powering source 
with a luminosity of 10® Lq= 3.9 x 10'*® erg s“*. Taking 
a CSM luminosity contribution of even 5 x 10®® erg s“*(see 
Section 1531) and an optical depth of « Ar « 2.3 yields 
a total Ha luminosity arising from CSM interaction of 
~ 5 X 10®® erg s“*. While this emission alone is not suf¬ 
ficient to heat the dust to the observed temperature, the 
bulk of emission from CSM interaction most likely arises in 
the X-rays and UV, where there are many emission lines 
at wavelengths where dust grains absorb very efficiently. 
Such an effect is common ly observed in supernova remnants 
(e.g., Temim et al.| 201 and other interacting SN e (e.g. , 
iFransson et al.ll2002l l. In fact, IChevalier &: Franssonl (1 199^ 1 
show that the conversion efficiency from X-rays to Ha is only 
~ 1%, suggesting X-ray luminosities of ~ 5 x 10** erg s“*. 

While we do not have any X-ray observations available 
at the time of this study, we do have nearly contemporaneous 
UV observations (see Table [T|). Taking into account an opti¬ 
cal depth of Td « Ar « 2.3 yields upper limits on the near- 
UV and far-UV fluxes of 5.5 x 10*® and 6.8 x 10*® erg s“* , 
respectively, for a total UV output of ~ 2.3 x 10** erg s“*. 

While the UV output is still not sufficient to power the 
K'-h&nd emission at the vaporisation radius, we note sev¬ 
eral caveats. Again, we stress that a significant portion of 
the flux arising from CSM interaction is likely emitted in 
the X-rays, for which we do not have data. Furthermore, 
the vaporisation radius sets only a lower limit on the dust- 
shell radius. Figure [T^j shows that a shell with a radius just 
a factor of 3 larger can yield the same iV'-band flux with 
a more massive dust shell that is > 100 K colder, which 
in turn requires ~ 10® Lq= 3.9 x 10** erg s“* to heat it. 
This scenario is consistent with the upper limits provided 
by the UV observations and X-ray estimates (assuming a 
1 % conversion efficiency). 



Figure 14. Velocity profile of the Ha line. Although noisy, the 
Ho line is broad (FWHM = 2000 ± 200 km s“^). At the very 
least, this rules out the possibility of narrow emissio n from an 
unde rlying H H region dominating the line flux le.g.. I^x et al.l 
l2013h . A light-red bar signifies pixels flagged for having bad data 
quality. 


6 DISCUSSION AND CONCLUSION 


In this article, we present new observations of the 3000 day 
plateau in SN 2006gy, including Keck/NIRC2-AO-LGS K'- 
band photometry on day 3024 and an HST/STIS spectrum 
on day 3001. After examining the energetics, colour evolu¬ 
tion, dust mass, and presence of broad Ha in the optical 
spectrum, we find that the new data are consistent with 
the presence of both a scattered optical light echo and IR 
emission from dust radiatively heated by X-rays and UV 
emission from ongoing CSM interaction. 

Even considering the 3000 day plateau, the total radi¬ 
ated energy output calculated in Section [5.21 is only Etot = 
2.9 X 10®* erg, which is not significantly greater than the 
total energy emitted by the SN during the first ~220 days 
an d still within t he 5 x 10®* erg of kinetic energy calculated 
by ISmith et al.l (120101 1. The derived energy budget there¬ 
fore does not contradict the CSM interaction model for SN 
2006gy that inv okes a large CSM formed by a very massive 
progenitor star ([Smith et al.ll2010fl . 

Assuming that the CSM interaction scenario provides 
the dominant power source for the day 3000 IR dust emis¬ 
sion, the total dust-shell mass can be tied to the progenitor’s 
total mass-loss rate, 


M, 
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M© yr 


(5) 
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for a dust-to-gas mass ratio = Ma/Mg « 0.01, and a 
progenitor wind speed Hw The relatively narrow lines ob¬ 
served in SNe Iln originate in the slow pre-shocked CSM 
and can be used to approximate the progenitor wind speed. 
The narrow lines observed in SN 2006gy have w idths corre¬ 
sponding to Ww = 200 km s“^ JSmith et al.ll200^ . Assuming 
a thin shell Ar/r = 1/10, dust-shell radius r = 3 x 10^^ cm 
(see Section FS. II) . and dust mass Ma = 0.1 Mq, the approx¬ 
imate mass-loss rate to produce the observed dust shell is 
A/outer ~ 0.2 M© yr“^. A larger shell radius would require 
a larger dust mass (see Figure [T2|, but these nearly cancel 
each other out and therefore do not change the calculated 
mass-loss rate significantly. 

The optical and/or UV/X-ray emission generated by 
CSM interaction can served as a tracer for the mass loss at 
the inner, shock radius. Assuming a density oc r~^ wind pro¬ 
file, the rate can be written as a function of the optical/X-ray 
luminosity, progenitor wind speed, and shoc k velocity (e.g., 
IChugai fc Danzig^ll994l : ISmith et al.ll2009l '): 


Mn 


2uw - 

-^Topt/UV/X, 

2.1 X 10 

-1 
e 

0/5 


^opt/UV/X 


3 X 10"^! erg s 

Uw 


120 km s"i 
-3 


104 km s 


-1 


M© yr 


( 6 ) 


where e < 1 is the efficiency of converting shock kinetic 
energy into light. We assume a value e « 0.1, although 
the conversion efficiency can vary with wind density and 
shock speed. Depending on the extinction correction we ap¬ 
ply, a late-time optical luminosity of Lr ~ lO^^ erg s“^, 
wind sp eed Uw = 200 km s”^^ shock velocity Vs = 4000 
km s“^ l|Smith et al.ll200^ . and conversion efficiency e = 0.1 
correspond to a mass-loss rate Minner ~ 10“^ M© yr“^. 

Equations [5] and |6] yield only order-of-magnitude ap¬ 
proximations, but it appears that the inner CSM mass-loss 
rate is an order of magnitude smaller than the mass-loss rate 
for the outer dust shell. The dust shell was likely formed dur¬ 
ing a period of increased, nonsteady mass loss. Furthermore, 
assuming spherical symmetry, the total implied dust mass 
from Figure \T2\ is > 8 x 10“^ M©. While this number is 
comp arable to that of many other SNe Iln iIFox et al.ll201ll . 
l2013h , it represents some of the most massive shells observed 
around SNe Iln. 

As a comparison, consider that the total warm dust 
mass of SN 2006gy is nearly as large as the total cold 
(~ 20 K) dust mass observed in the extremely dusty SN 
1987A, whic h was detected only due to its nearby distance 
from Earth llMa tsuura et al.l 1201 ll : lindebetouw et al] l2014l : 
iMatsuura et al.l 2015l l. The implication is that much larger 
dust reservoirs may be hiding in SN 2006gy. If this much dust 
forms in other SLSNe and can survive the forward and/or re¬ 
verse shocks , it may be able to account for the cosmic dust 
budget (see iGall et al] l2014l . and references therein). The 
large dust shell and extremely dense CSM may also explain 
the suppressed X-ray and radio emissi on, an effect observed 
to lesser degree s in other SNe Iln (e.g., IVan Dvk et al.lfl9^ : 
IFox et ^l2009l l. 

The total dust mass of > 10“ ^ M© derived from CSM 


interaction indicates a total shell mass of ~10 M©, assuming 
stand ard g as-to-dust mass rat ios. As noted bv ISmith et al.l 
ll2008ll and I Miller et al.l (l2010l ~) , such a large envelope sug¬ 
gests a very massive progenitor star . These numbers ar e sur- 
prisingly clo s e to t hose derived by ISmith et al.l (l2008l l and 
I Miller et al.l (l2010ll . which is probably not a mere coinci¬ 
dence. More likely, the same du st sh ell observed as a th er- 
mal echo by ISmith et al.l ll2008ll and iMiller et al.l ll2010h is 
increasingly heated by CSM interaction. 

The derived radius of the dust shell in this paper, how¬ 
ever, d iffers from that of ISmith et al.l ll2008h and llVIiller et al.l 
ll201flll . For the CSM interaction scenario, the shell resides 
nearly a factor of 10 closer, which would suggest the ob¬ 
served dust shell was formed in an outburst that occurred 
as recently as just a few hundred years prior to the explo¬ 
sion. The time frame of the outbursts may offer important 
constraints on the progenitor mass-loss mechanisms. As the 
forward shock overtakes this shell of material, we expect the 
SN luminosity to increase for a period of time consistent with 
the width of the shell of material. Continued monitoring of 
SN 2006gy will help to constrain the dust-shell properties. 

Overall, this paper highlights the significance of multi¬ 
wavelength observations, even if limited in wavelength range 
and sensitivity. Future observations of interacting SNe 
should cover more wavelengths (e.g.. X-ray, mid-IR, and ra¬ 
dio) over more epochs. Ongoing observations of SN 2006gy 
will ultimately constrain the precise extent of this material, 
the presence of additional shells, and the possibility of al¬ 
ternative, more exotic energy sources, such as a magnetar. 
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